Dipolar effect in coherent spin mixing of two atoms in a single optical lattice site 
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We show that atomic dipolar effects are detectable in the system that recently demonstrated 
two-atom coherent spin dynamics within individual lattice sites of a Mott state. Based on a two- 
state approximation for the two-atom internal states and relying on a variational approach, we have 
estimated the spin dipolar effect. Despite the absolute weakness of the dipole-dipole interaction, it 
is shown that it leads to experimentally observable effects in the spin mixing dynamics. 
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The recent successes of coherent spin mixing Q, 1^ 
inside Bose condensed atoms, driven by reversible colli- 
sions between pairs of atoms |Afi? = 0) + \Mp = 0) ^ 
I — 1) -I- 1 1), have generated significant interest in the 
quantum dynamics of atomic spins 4, 5,, ^, These 
experiments cover a broad limit from a condensate with 
large number of atoms P, to many identical trapping 
sites as in an optical lattice each containing two atoms 
0, and raise significant hope for the long discussed ap- 
plications of atomic quantum gases to the emerging field 
of quantum information science. As observed in these 
experiments, the spin coherence time rivals the best mo- 
tional state coherence ever achieved on neutral atoms Q 
and is ideally suited for quantum information processing 
applications I9Hiq|. 

In this paper, we reveal an interesting observation: 
the atomic spin dipolar effects seem to be detectable 
in the system that recently demonstrated two-atom co- 
herent spin mixing Dipolar interactions are ubiq- 
uitous in atomic systems. It offers a unique play- 
ing field because of the two significant differences from 
the the nominal s-wave short-ranged interactions: it is 
anisotropic and long-ranged [Tlj . These interesting prop- 
erties have stimulated intensive research of dipolar con- 
densates mill El [mil. 

Our model system consists of two spinor atoms (^''Rb) 
inside a harmonic trap. Corrections due to the anhar- 
monic trapping potential of a standing wave optical lat- 
tice E3 will be addressed elsewhere. We will consider 
the simple case of a cylindrical harmonic trap instead of 
the spherical harmonic trap as in the experiment Q , be- 
cause the dipolar interaction is known to display a sen- 
sitive dependence on the trap aspect ratios The 
spherical harmonic trap minimizes the dipolar effect and 
thus provides an excellent calibration for the experimen- 
tal systems, particularly the small spin exchange interac- 
tion n], fl?! . Following the usual procedure of separating 
the total motion into the center of mass (CM) and the 
relative (rel) motion, the system Hamiltonian becomes 
H = i/cM + ^^rci with 
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for the CM- and rel-motion, respectively, with coordi- 
nates R = {ri + f2)/2 and r — ri — r2- M = 2m is the 
total mass, while /i = m/2 is the reduced mass. 7?rei 
contains several parts as outlined below. 
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where Hq describes the rel-motional, and accounts for 
the regularized s-wave contact interactions between two 
spin-1 atoms with 6'^''^^\r) = 5{r){d / dr)r. Hdd denotes 
the spin dipolar interaction. The first order Zeeman ef- 
fect does not contribute due to a zero magnetization, 
while the second order Zeeman effect Hb ~ ±72 (F • BY 
(Hz/Gauss^) in reference to the Mp ~ state with the 
'-f ' and '— ' signs for F = 1 and F = 2 respectively. The 
various interaction coefficients are listed below 
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ao(2) is the scattering length for the combined channel of 
total F — 0(2). gp is the Lande g-factor for the hyperfine 
spin state of F, and /is is the Bohr magneton. 

Before presenting our theoretical analysis, we comment 
on the strength of spin dipolar interactions ^3. In a 
scalar condensate, the dipolar effect is usually calibrated 
against the nominal s-wave interaction, i.e., one simply 
compares \cd\ with |co| and uses their ratio as a param- 
eter. This ratio can be increased by decreasing cq per- 
haps through a Feshbach resonance ll9'|. For spin dipoles, 
however, the spin exchange interaction coefficient C2 also 
needs to be compared. The extremely small C2 for ^^Rb 
serves to enhance the spin dipolar effect because the 
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spinor nature of the ground state energy is determined 
by C2 not by the much larger cq. Using |cd|/|c2| as a pa- 
rameter, ^^Rb could be viewed as a stronger spin dipolar 
condensate than Cr. Although in a polarized conden- 
sate, Cr atoms {F = 3) with a dipole moment of 6fiB 
enhance significantly the dipolar interaction. When the 
spinor nature of the dipole is of interest [isl I20I I2H , the 
dipolar effect for Cr is weaker because of the large spin 
exchange interactions [isl |2^ . 

We first will estimate the relative strength of different 
interaction terms for the situation as in Ref. For 
an optical lattice oc Vo[sin^(fca:) -l-sin^(fcy) -I- A^sin^(A;z)] 
of a depth Vb = sEr in units of the recoil energy Er 
for the lattice laser, each single lattice site is approx- 
imated like a harmonic trap with a radial frequency 
ujp = y^2sErk^ /m = 2y/sEr/h and a Gaussian relative 
motional ground state 
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with a radial width Up — yJJtJJjjIjjp) {p = ^yx^ + y^). 
The center of mass motional ground state ^q{R) takes 
the same form except for the width Ap = h / [2mLU p) . 
For a wavelength of 840 nm, we find ujp ~ (27r)41, 126.3 
Hz, and Op = 0.0752 /im for s — 40. We further note that 
ao = 101. Sob and 02 = 100. 4aB, with as being the Bohr 
radius These give rise to typical density interac- 

tion co{n)o/h ~ (27r)6, 589 Hz, spin exchange interaction 
C2{n)()/h ^ (27r)(— 30) Hz, and spin dipolar interaction 
Cd{n)o/h ^ (27r)2.71 Hz for the spherically symmetric 
case of A = 1, all much less that the trap level spacing. 
Thus, a reasonable estimate of the motional state would 
be the non- interacting ground state 4>o (r) from which the 
averaged density becomes 
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We further note the relative strength of C2/C0 ~ 
—0.00462657 and Cd/c2 ~ —0.0902, which encourages a 
Gaussian variational ansatz ^12j . 

The dominant mixing interaction is spin exchange that 
couples the two-atom internal state \Mp — 0, M'p = 0) 
to 1 1,-1). The spin dipolar interaction is about 5 
times smaller and averages to vanishingly small net ef- 
fect for a spherical symmetric motional state. We there- 
fore will limit our discussions to the above two internal 
states, a picture uniformly adopted by the experimen- 
talists |li ^, 3J. We have further carried out numerical 
simulations at _B = of the full system dynamics, includ- 
ing all other spin states that are coupled by the dipolar 
term, i.e., with the complete Hilbert space of spin de- 
gree of freedom: \Mf = -l,0,l;Af^ = -1,0,1). For a 
spherical trap with the same (averaged) frequency, the 
probability of atoms in other spin states that do not con- 
serve the total magnetization (|1, 1), | — 1,— 1), |0,l),and 
|0, — 1)), which are responsible for the dipolar relaxation. 



is found to be less than 10~^ in the first oscillation period 
of coherent spin mixing; even for a trap with A = 3, this 
probability remains negligible, and is in fact only sev- 
eral times enhanced. Thus the negligibly small dipolar 
spin relaxation in ^^Rb 23] makes the two-state model 
an excellent approximation for our system. At finite val- 
ues of the B-field, except for accidental resonances when 
other spin states in higher motional states are shifted into 
near resonance with the two-state doublet in ground mo- 
tional state, the linear Zeeman effect generally leads to 
large detunings, also validates the approximation. Even 
at accidental resonances, the total population out of the 
two-state doublet is found to be only ~ 10~^ for a spher- 
ical trap. For Cr atoms, however, more effort is needed 
to understand the conditions for spin mixing dynamics 
due to the much enhanced dipolar relaxation 24] . 
The two-atom wave function is then approximated as 

ao,o|0, 0)V'o,o(^i,n2) + ai^-ijl, -l)V'i -i(?^i, ^2)- (H) 
This leads to a spin mixing matrix element of 

hn = I df, J df2roAn,r2){0,0\iHs 

+ifdd)ll,-l)^i,-i(ri,f2). (12) 

As a first approximation, the relative motion is sim- 
ply taken to be the ground state of the harmonic trap, 
i.e., '0AfF,M;,(n,r2) = ^o{R)(t>MF,M'^ir) with (/)o,o(r) = 
(j)i.-i{r) = <j)a{r). This leads to 
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df^o o(^)(0, 0\H, + Hdd\l,-l)^i.~i{r). (13) 



An improved approximation is the variational calcu- 
lation, labelled as 4>^'"\ for the relative motional state 
including the dipolar interaction as have been used ex- 
tensively in the past We take a Gaussian ansatz 
with its widths Wp/z as variational parameters |3l 
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withx(K) = [2K^ + l-3K^H{K)]/[2{K'^-l)] + {fi'^-l)H{K) 

and H{k) = taii]i^^\/l — k^/\/1 — k^. = ^yh/(pLUz) 
is the axial width of the trap, and n = Wp/wz is the 
aspect ratio of the variational ground state, non-spherical 
(or K y I) even in a spherical harmonic trap with A = 
ijJzlijJp = 1 because of the dipolar interaction [12]. 0o,o('^) 
is obtained from the result of (^i^_i(r) by excluding the 
dipolar interaction or taking = and adjusting to its 
own s-wave scattering strength by taking C2 = 0. 
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A mean field approach is sometimes used in the Htera- 
tm'e where the two-atom motional state '4'Mp.M'^i^ii^2) 

is approximated by '^^ {fi)(f>c^'' {f2) as for a two- 
atom condensate with 0c^^'*^^(r) obtained from the cor- 
responding Gross-Pitaevskii equation 
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with Vint = Co\(j>c{r)\'^ , or 
Vint = (co-C2)|0c(r)|' 
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respectively, for the internal states |0,0) or |1,— 1). 
is the chemical potential, and 9 is the angle between the 
z-axis and r — ? . A Gaussian variational approach is 
proven to be adequate within the parameter regions of 
interest |l2l \v\ . I25II. Like the 0o approximation, atom- 
atom correlation [2g is neglected because of the use of 
product motional states here. 

We have used the momentum space pseudo- 
potential v{kS') = -yi' l{2Tx'^m)o!^^sh{P-i{c^s.ek>) + 
{k/k'fP2{cos9k)] with a'^^ = V2mcd/(12v^?l^) for a 
more accurate evaluation of the dipolar term in Eq. (|14|l . 
which can differ upto 50% for the parameter range re- 
ported here. Because the complete motional wave func- 
tion is Gaussian shaped, we find that {(l>i^-i\Hdd\(t>o.o) 
can be evaluated analytically. In the results shown 
below for two ^^Rb atoms, the spin mixing effective 
Rabi frequency is defined as ^IcS — vrP + A^, and 
fes = flcs/^TT. A is the bare energy diflference: A = 
{Ei^-i - Eofl)/h. 




log^o(^) 



FIG. 1: (Color online) The aspect ratio dependence of f^s 
computed with the three approximation schemes as labelled. 
The three solid lines include dipolar interaction, while the 
dot-dashed lines are the corresponding ones without dipolar 
interaction. Here we fix [lo^lu^Y^^ = (27r)41.1(kHz). 

In Fig. m we show the dependence of /eft on the trap 
aspect ratio A at an external magnetic field of B = 0.1 



(Gauss). To facilitate a fair comparison, we have fixed 
the geometric average of the trap frequencies (w^^)^/'^ to 
the spherical trap oiujp = = (27r)41.1(kHz) 0. Some- 
what surprising at first sight is the noticeable quantita- 
tive differences (within experimental sensitivity) among 
the different approximations. At A = 1 for a spherical 
trap, our results also differ from the experimental obser- 
vations We note that the harmonic approximation 
of Vbsin^(A:a;) by muj'j,x'^/2 introduces about a 5% error 
[l7|. Although the spread among the different approxi- 
mation schemes calls for a more accurate treatment for 
the motional state, the dipolar effect due to the Hdd term 
in Eq. IjlSI) is quite accurately reproduced to less than 
1% 28]. An improved treatment of the relative motional 
wave function is also needed inside a cylindrical trap, 
if spin mixing is used to calibrate atomic interactions 
like what has been accomplished for a spherical harmonic 
trap [13 ■ 
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FIG. 2: (Color online) The B-field dependence of /off com- 
puted within (j/)*^"' approximation for A = 3. Solid line includes 
dipolar interaction while the dot-dashed line excludes dipolar 
interaction. Here, unlike in Fig. we fix ujp — (27r)41.1(kHz). 

Based on the current experimental sensitivity, dipolar 
effects should be detectable for A > 3 and can constitute 
a more than 10% increase in fgg. They are minimized for 
a spherical trap as shown in Fig. |21with an actual shift 
of about 10"'' or less for the experiment of Ref. [T^ . 

Also we have studied the F = 2 manifold where 
spin mixing dynamics is generally dominated by three 
two-atom internal states |0, 0), |1, — 1), and |2, — 2) at zero 
magnetization We find that spin dipolar effect re- 

mains observable in the frequencies of the various mixing 
channels. More interesting is the two-state mixing chan- 
nel of I — 1, —1) ^ |0, —2) with a nonzero magnetization 
as shown in Fig. 13 Inside an oblate trap, dipolar interac- 
tion is dominantly repulsive in state | — 1, — 1) in contrast 
to the attractive state |1, —1) of the F = 1 case. At a 
weak magnetic field and for A = 3, we find the dipolar 
interaction constitutes a 25 Hz downward shift computed 
within the 0^"^ approximation. 

Dipolar interaction is normally stronger by a factor of 
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FIG. 3: (Color online) The same as in Fig. |5| except for spin 
mixing between | — 1, —1) and |0, —2) for F — 2. 

two in magnitude along the direction of the dipoles (fa- 
vored by the A < 1 geometry) in comparison to the per- 
pendicular direction (favored when A > 1). Our results 
on the spin mixing frequency, however, reveals a com- 
pletely opposite trend; we find relatively larger (smaller) 
dipolar effects for A > 1 (< 1). This can be easily under- 
stood. For the F — 1 case, the dipolar term Hdd is pos- 
itive (negative) for A < 1 (> 1), thus destructively (con- 
structively) add to the (negative) C2 term in Eq. (|13|l . 



For the F = 2 case, the same reasoning applies despite of 
the opposite dipolar interaction in state | — 1,-1). The 
dipolar effect becomes constructively enhanced for A > 1 
because the C2 is positive in F = 2. 

Before conclusion, we note that we have also calcu- 
lated the spin mixing frequency for two ^^Cr atoms with 
F = 3. In the limit of very large (A — 10) and very small 
aspect ratios (A = 1/10), dipolar interaction causes about 
an 8% difference in the spin mixing frequency, which is 
about four times smaller than in F = 1 of ^^Rb at A = 10, 
and close to each other at A = 1/10. Similar conclusions 
hold for ^^Rb in F = 2. It is in this sense we say that the 
dipolar effect of ^^Cr is weaker than ^^Rb. The absolute 
frequency difference due to dipolar interaction, however, 
is larger in ^^Cr because of the faster spin dynamics in- 
volved due to the large spin exchange term in ^^Cr. 

In summary, we have studied dipolar effects in spin 
mixing between two atoms trapped in a single optical 
lattice site. While this effect is small, and can be ignored 
completely for spherical harmonic traps, we find it is ob- 
servable inside cylindrical harmonic traps, especially for 
oblate shaped traps with A > 3. We hope this study will 
stimulate experimental efforts aimed at observing dipolar 
effect in spin mixing. 
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